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Beginning with an effective field theory based upon meson exchange, the 
Bethe-Salpeter equation for the three-particle propagator (six-point function) 
is derived. Using the one-boson-exchange form of the kernel, this equa- 
tion is then analyzed using time-ordered perturbation theory, and a three- 
dimensional equation for the propagator is developed. The propagator con- 
sists of a pre-factor in which the relative energies are fixed by the initial state 
of the particles, an intermediate part in which only global propagation of the 
particles occurs, and a post-factor in which relative energies are fixed by the 
final state of the particles. The pre- and post-factors are necessary in order to 
account for the transition from states where particles are off their mass shell 
to states described by the global propagator with all of the particle energies 
on shell. The pole structure of the intermediate part of the propagator is used 
to determine the equation for the three-body bound state: a Schrodinger-like 
relativistic equation with a single, global Green's function. The role of the 
pre- and post-factors in the relativistic dynamics is to incorporate the poles 
of the breakup channels in the initial and final states. The derivation of this 
equation by integrating over the relative times rather than via a constraint 
on relative momenta allows the inclusion of retardation and dynamical boost 
corrections without introducing unphysical singularities. 

21.30.Fe, 21.45,+v, 13.75.Cs 



I. INTRODUCTION 



Faddeev's three-body work of 1960 JTJ, combined with the earlier two-body work of Bethe 
and Salpeter, and Gell-Mann and Low produced significant interest within the nuclear 
and particle physics communities in solving the fully relativistic three-body problem 



Calculations based upon the full four- dimensional theory, however, have proved very difficult, 
and have only recently been performed by Rupp and Tjon using separable interactions |TT 
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In the meantime, three-dimensional nonrelativistic calculations based upon the Schro- 
dinger equation have progressed significantly. Using NN potentials which provide a good 
description of two-body data, recent calculations have been performed to a precision of 
10 keV ||12|| . Given the successes of these nonrelativistic calculations in obtaining excellent 
precision, discrepancies with experiment of order 10 keV or larger are due to inaccuracies 
in the theoretical input rather than uncertainties in the calculations. The main "missing 
physics" in these calculations are three-body forces and relativistic effects. A recent calcula- 



tion found the triton to be underbound by roughly 480-860 keV [12]. A recent assessment of 



the triton potential energy suggested that a consistent relativistic calculation could account 
for as much as 300 keV of repulsion JO}, and references therein]; however, recent relativistic 
calculations using the Blankenbecler-Sugar formalism and the CD-Bonn potential showed 



an increased binding of 200 keV |i~2" |. The size of the relativistic effects determines the 
amount of three-body interaction that would therefore be required. A consistent relativistic 
three-nucleon calculation is needed in order to permit some understanding of the respective 
roles of relativistic effects and three-body forces in nuclear binding. 

In view of the successes of the nonrelativistic calculations and the need for relativistic 
calculations, many three-dimensional reductions of the relativistic four- dimensional Bethe- 
Salpeter equation have been made following the two-body work of Blankenbecler and Sugar, 
and Logunov and Tavkhelidze [14-19]. This approach, often called the quasipotential ap- 
proach, consists of replacing the Bethe-Salpeter equation with a set of two coupled equations. 
These equations involve two new functions, the quasipotential W and the quasipotential 
Green's function Gqp. We may choose one of these functions arbitrarily. The requirement 
that these equations be equivalent to the Bethe-Salpeter equation then fixes the other func- 
tion. Traditionally, the quasipotential Green's function Gqp is chosen to contain a Dirac 
5 function constraint which reduces the dimensionality of one of the new equations. Diffi- 
culties with this procedure arise when the equation for the quasipotential W is truncated. 
What form these difficulties take is dependent upon the form of the constraint used in the 



quasipotential Green's function. For example, in the Gross equation |2(J this truncation 
introduces into the wave function unphysical singularities which must be removed by hand. 
In the instant formalism ||21|| , singularities arise when attempting the dynamical boost of 
the wave function. Although they take different forms, the root cause of these difficulties 
is the 5 function constraint combined with the truncation of the quasipotential. Given 
the difficulties with solving the full four- dimensional equation, the success of the nonrela- 
tivistic three-dimensional calculations, and the fundamental problems in the quasipotential 
approach, a different technique for dimensional reduction seems warranted. 

In this paper we develop an approach to the relativistic three-body problem, with an 
emphasis on three-body bound states. We start from the Bethe-Salpeter equation for the 
full three-body Green's function in momentum space. For definiteness, we use scalar me- 
son exchange as a model for the interaction; extension to other forms of interactions is 
straightforward. Negative energy states are omitted as the dominant physics is obtained 
from positive energy states. 

We perform a Fourier transformation of the zeroth component of all internal momenta 
to relative-time variables, and carry out the relative-time integrations, which has the effect 
of transforming each Feynman graph into several time-ordered graphs. These time-ordered 
graphs are three-dimensional in nature. Rearranging and summing graphs produces an 
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FIG. 1. The sum of two-body- irreducible graphs denoted by Vj. 



expression for the full, four-dimensional, three-body Green's function in which all of the 
internal variables are three dimensional. It follows that the bound-state equation is three 
dimensional as well. We find that the propagator consists of a pre-factor in which the 
relative energies are fixed by the initial state of the particles, an intermediate part in which 
only global propagation of the particles occurs, and a post-factor in which relative energies 
are fixed by the final state of the particles. The pre- and post-factors are necessary in 
order to account for the transition from states where particles are off their mass shell to 
states described by the global propagator with all of the particle energies on shell. This 
formalism allows calculations of bound states in three-dimensions (where much success has 
been shown) and provides the formalism for embedding the result within a four- dimensional 
covariant scattering theory. 

In Sec. |II|we define the full three-body Green's function. Then, in Sec. HJwe examine the 
three-dimensional reduction of the internal momenta of the Green's function. In Sec. [IV] we 
examine the structure of the Green's function, organizing the summation of graphs into pre- 
and post-factors, and a three-dimensional iterative Green's function. In Sec. [V] we extract 
the bound-state equation from the pole of the Green's function. Finally, in Sec. [VI] we discuss 
our conclusions from this work. We also provide three appendices: Appendix [A| in which we 
provide the rules for time-ordered perturbation theory for our model, Appendix [B], in which 
we provide more details of the reduction from four to three dimensions, and Appendix 0, 
in which we discuss cluster separability in our formalism. 



II. THE FULL THREE-BODY GREEN'S FUNCTION 

The four-dimensional three-body Green's function is defined in field theory as the six- 
point function 

G(x 1 ,x 2 ,x 3 ;y 1 ,y2,y 3 ) = (0| T[ip 1 (x 1 )^ 2 (x 2 )^3(x 3 )^ 3 (y 3 )^2(y2)'ipi(yi)) |0) . (1) 

Here we allow ip to represent either spin- 1/2 or spin-0 particles, and consider the three 
particles to be distinguishable. Note that with distinguishable particles, ipi can only contract 
with ipi, and for spin- 1/2 particles anticommutes with ip2 and ip 3 . The interaction is assumed 
to be a sum of meson-nucleon interaction terms [see Eq. (g)]. 

The Bethe-Salpeter equation |2||| for the three-body Green's function may be derived 
from Eq. ([!]) by expanding Q in a perturbation series in the interaction picture. Rearranging 
the resulting Feynman graphs into two sets, the two-body irreducible graphs and the iterative 
graphs, and denoting the sum of two-body irreducible graphs for cluster j as Vj (see Fig. p, 
we find 

3 

G = -igi9293 Vj9k9i = Go + GJ2 [Vffai) -1 £o 

i=l j 

= Go + GoE[^93r l ]G, (2) 
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where Qq = —igig^gs- (For reviews of the Bethe-Salpeter equation see Nakanishi J22| and 



Remiddi ||23|| .) Note that we are considering distinguishable particles, and that the subscript 
of the single-particle Green's function gj refers to the particle species. For indistinguishable 
particles the products would need to be symmetrized or antisymmetrized depending upon 
whether we are considering spin-0 or spin-1/2 particles. The index j of the interaction in 
Eq. @ conforms to "odd particle out" notation, such that Vj represents the sum of two- 
body irreducible graphs between particles k and I, where (jkl) is an even permutation of 
(123). The "self-energy-summed" single-particle propagator gj is approximated by the free 
propagator using the physical mass, 

g^x'^^-mT^x')^)]^) . (3) 

Although this analysis neglects Feynman diagrams representing three-body interactions, we 
shall see later that Eq. (0) contains time- ordered diagrams representing three-body interac- 
tions. 



III. THREE-DIMENSIONAL REDUCTION 

A three-dimensional reduction follows from decomposing graphs into sums of distinct 
time intervals between consecutive meson emission and absorption events, and then inte- 
grating over the duration of each interval. 

In order to establish some notation, the free-particle propagator is written as the sum of 
a positive- and a negative-energy part 

M = (4) 
where the on-shell energy of a particle is denoted by 



and we define 



and similarly 



e,(p) = vV + m? (5) 



N + (p) = u(p)u(p), (6) 



N'(p) = -T,(-p)tj(-p) . (7) 

For spin-1/2 particles, u(p) is a positive-energy Dirac spinor, v(— p) is a negative-energy 
Dirac spinor, u(p) = w^(p)7°, and v{— p) = i>'(— p)7°. Dirac spinors obey the Hermitian 
normalization conditions: u'(p)u(p) = 1, and v*(—p)v(—p) = 1. For spin-0 particles 
w(p) = u(p) = l/v2e, and v(— p) = — v(— p) = \/y/2e. Combining these definitions with 
that of T provided in Eq. (PH|), we have Tu(p)u(p) = gom/e and Tv(— p)v(— p) = —g^mje 
for either spin-0 or spin-1/2 particles. This notation permits the analysis to proceed on a 
common footing for both spins. 

To facilitate the time-ordered analysis, we perform a Fourier transformation of the time- 
like component of the momentum of g(p), 
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= (-i)0(f - t)N + (p)e- ie{ ^ t '- t) - i6(t - t / )iV"(p)e' je(p){ *'-* ) . (8) 



In order to simplify the analysis, two approximations are used: 

(i) Positive energy particles. The contribution of negative-energy states to the propagator 
of Eq. (H) is neglected. This is the traditional starting point in nuclear physics. It breaks 
the covariance of the theory but is generally believed to be commensurate with our present 
understanding of nuclear forces. Extension of the analysis to incorporate negative-energy 
states will be left to the future. 

(ii) One boson exchange. For NN interactions, a suitable effective interaction should 
describe the NN phase shifts and deuteron binding. This may be accomplished by the use 
of a one-boson-exchange (OBE) interaction, such that Vj is replaced by the first term in 
Fig. [p. 

Consider an interaction Hamiltonian of the form, 



Hi 



~9o ■ ipjCftpj ■, for spin-1/2, 



j t (9) 
—2i7ijg : Q*ip&j :, for spin-0. 

I 3=1 



where go is a coupling constant. Separate fields are introduced for each particle in order to 
treat them as distinct particles. 

It is convenient to define a vertex factor as follows: 

U, for spin-1/2 
J y ^rnjgo, tor spm-U, 

where the factor 2mj is introduced in order that both cases have a common nonrelativistic 
limit. Performing a Fourier transform with respect to the timelike component of the mo- 
mentum, and using the approximation in which only the one-boson-exchange potential is 
retained, produces 

(-i)9(t k - t t ) — — — + (-i)0(t, - h) — — — | , (11) 

where u = \/// 2 + q 2 , with /i being the meson mass and q being its three momentum. 

The OBE interaction together with the use of the free-particle propagator [Eq. @] means 
that the analysis pertains to the ladder approximation Feynman graphs . 

Inserting Eqs. (|8|) and (|ll|) into Eq. (|2|) leads to an equation for the Green's function in 
which all of the "internal" integrations are over the three momenta and the relative times. 
Each contribution can be decomposed into a sequence of time intervals between consecutive 
meson emission and absorption events. The integral over the time duration of each such 
interval may be performed analytically to produce the time-ordered perturbation theory 
(TOPT) rules and corresponding graphs (see Appendix [A]). 

In order to provide a few examples, consider subgraphs that contain no initial or final 
particle lines, i.e., that are embedded within other graphs such that interactions separate 
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them from initial or final particle lines. Time intervals in which three particles propagate 
freely correspond to 



P u - ei - e 2 - e 3 + ir] 

and time intervals beginning with emission of a meson and ending with its absorption cor- 
respond to 



r 2 iV 2 + 1 TxNf 

V2u3 -P - ei - e 2 - e 3 - o; 3 + ir] ^ET 3 



(13) 



When a second meson is in flight during a time interval such as in Eq. flT3"|), the contribution 
is modified to the form, 



oc 



_ t 2 n£ 1 TxN? 

~— ____ v / 2^3 -P - ei - e 2 - e 3 - u x - uj 3 + irj y/2uJ^ 

(14) 

where labels of the exchanged bosons correspond to "odd-particle-out" notation. In Eq. (|14|), 
the expression on the right-hand side corresponds to the time interval and vertices within 
the dashed box on the left-hand side. For each distinct time interval there is a denominator 
equal to the total energy minus the sum of the on-shell energies of all of the particles present 
during that time interval. Although consideration was limited to the OBE interaction, the 
time-ordered rules apply quite generally and, for example, the "cross-box" diagram shown 
yields: 



rii\tf, i 



2uj' 3 p0 - e " - £ 2 - e 3 - u' 3 + ir) 



Y 2 N± 1 T 2 N+, 



P° - e'{ - e'i - e 3 - uj 3 - u' 3 + it] 

i ryvj 

P° - e'{ - e' 2 -e 3 -LU 3 + ir] ^ET 3 



(15) 



Note that there is an implied integration over the loop three momenta. In general there are 
six time-ordered diagrams corresponding to a single cross-box Feynman graph, differing by 
the time ordering of vertices on one particle line with respect to those on the other particle 
line. 

The transition from four to three dimensions necessitates a reclassification of diagrams 
in terms of three-particle irreducibility with respect to G so as to distinguish between 
two- versus three-body forces in the time-ordered formalism. Consider the simple, iterative 
Feynman (four-dimensional) diagram shown at the left of Fig. 0. The reduction from four 
to three dimensions produces 12 time-ordered diagrams, the sum of the first four being the 
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A a B b A b B 



+ 



FIG. 2. By the reduction from four to three dimensions, the Feynman diagram on the left 
produces the 12 time-ordered diagrams on the right. Note that four of the time-ordered diagrams 
are iterative like the original Feynman diagram, but that the other eight are not, and constitute a 
three-dimensional three-body force. 

iteration of the two-body force. The other eight, however, do not have the requisite time 
interval with only the three particles propagating, and are hence three-particle irreducible. 
These diagrams, along with their "sibling" diagrams (where particle two interacts with 
particle three before particle one), make up the three-body force Vg 31 . The subscript denotes 
"zero particles out" while the superscript labels the odd-particle-out exchanges of which this 
graph is composed. Please see Appendix |B| for a fuller description of this transition from 
four to three dimensions. 



IV. THE THREE-DIMENSIONAL THREE-BODY GREEN'S FUNCTION 

Having carried out the reduction to three dimensions and reclassified the diagrams, we 
may re-sum the infinite series for Q. This sum may then be separated and factored. First, 
we take cluster separability into account, and then we factor the connected three-body term 
into three pieces. 

Cluster separability [pi, 25(1 means that a Green's function describing the propagation of 



clusters of particles, when there are no interactions between the clusters, consists of a product 
of independent factors, one for each cluster. Each cluster's factor is the same as if the other 
clusters were not present. For three particles, we may have (a) no particles interacting, (b) 
two particles interacting and the third a spectator, and (c) all three particles interacting. 
As these are the only possible cases, cluster separability states that for the most general 
Green's function we have 



7 



G(pif,P2f,P3f,Pli, P2i , P&) = S^i (pif, Pli) Q% (P2 f\Pli) Q% (P3 / ; P3i 



3=1 

+ Gi 3) (Plf,P2f,P3f,PlhP2hP3i) • (16) 

We are denoting a fully connected n-particle Green's function for particles j, . . . , k by (?cy\.. 
(We omit the 1,2,3 label for n = 3.) Note that the fully connected part, Gf?\ vanishes by 
definition in the limit that one particle does not interact; the correct cluster limit obtains 
from the disconnected diagrams in this case. Note also that if a three-body bound state 
exists, it corresponds to a pole in Qj® and has no contribution in the other parts. Please 
see Appendix for a more complete discussion. 

As we are primarily interested in obtaining the bound state, we discard two sets of 
graphs: (a) disconnected graphs in which one or more of the particles never interact, and 
(b) graphs in which there is an initial or final particle in every three-body-reducible time 
interval (i.e., in which there are no factors of Go)- (See Appendix [C].) After discarding these 
two sets of graphs, we find that the sum of connected graphs in which there is at least one 
fully internal time interval is factored into three parts, 

— > i(2*)W\P f - PdfSr 3B) G { s B) f^ m) ■ (17) 

We shall define G%, / pre , and / post shortly. The superscripts (4D — > 3D), (3D), and (3D — ► 
4D) label the dependence of the three parts of the Green's function on the initial-state and 
final-state relative energies. The superscript (3D) denotes that G 3 has no dependence upon 
the initial or final relative energies. As we have integrated out the internal relative energies 
this function is purely three dimensional. The superscript (4D — ► 3D) on / post denotes 
that it depends upon the final relative energies of the particles. This function describes the 
transition from the three-dimensional internal time intervals to the four dimensional final 
state. The function / pre performs a similar function, depending upon the initial relative 
energies of the particles. As mentioned in Appendix |C|, this pattern of energy dependence is 
present in any fully connected Green's function that goes beyond the Born terms. However, 
the functional forms of the parts depend upon the number of particles involved. 

We define G3 as 



= Go 



ra=0 



E YiN^N? + V N+N+N+ ] Go 
\j=i 



= G + Go \p V3N+N+ + VoN+N+N+j G 3 . (18) 

The three-body potential (V& = V 12 + V 23 + V 31 + V 123 + • • •) is the sum of diagrams in 
which (a) all three particles are interacting and (b) there are no three-particle-reducible time 
intervals. The three-dimensional three-body Green's function G3 represents the sum of the 
diagrams in which all of the particles have interacted at least once, and will interact at least 
once more. This requirement ensures that interactions in G3 are separated by the global 
propagator Gq and do not contain any dependence on the initial-state and final-state particle 
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relative energies, i.e., are fully three dimensional. When there is a three-body bound-state 
pole, it is contained in G 3 . 

Additional factors / pre and / post arise in Eq. ( |iT|) from sequences of two-body interactions 
of particles j and k in the initial and final states, separated by the two-body propagator for 
cluster I, G l = 1/ [P° k - 6j - e k + ir] 
These factors may be expressed as 



where P° fe = p° + jP k is the total energy of the pair. 
NfN+N+ 



post 



(pp 



[ f - ei + iv)ip% ~ e 2 + iv)(p%f - e 3 + ^) 

' V 



j=i 



(19a) 



pre 



i 

S (Pi< - e i + «7)(p§i - e 2 + iv)(P°3i - e 3 + ^) ' 



(19b) 



where ilf and fi^ are the left and right two-body wave operators for cluster j, defined as 



m=0 

OO 



Pfc+Pi 
1 



m=0 



e k - ei+ irj 



e k ~ e t + if] 



E {ViNtNfGi 



m=0 



(20a) 
(20b) 



These wave operators transform the two-body free propagator Gq into the full two-body 

r 2- 



Green's function Go: 



G*2 — GqVL^ — VL^Gq 

= Gi + GiV j N+N+Gi. 



(21) 



Through these wave operators, the pre- and post-factors (/ prc and / pos t) contain the two- 
body bound-state poles for the different possible clusterings of particles within the initial 
and final states. Note that when all particles in the initial and final states are on the mass 
shell, Pj k = P° — e/, and G l = G . For consideration of interactions of the bound state with, 
say, a photon, one needs the full structure of Q, including the pre- and post-factors that 
allow for breakup. However, the three-body bound state is determined from consideration 
of G 3 as defined by Eq. (0). 



V. BOUND-STATE EQUATION 

Assume that there is a pole in G 3 at P° = E B (P), where E B (P) = \J + P 2 , and M B 
is the bound-state mass. To find the bound-state equation, we write G3 as 

G 3= WW, + R, (22) 

P u — E B + ITj 
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and therefore 



^ 3 )=z(27t)V 4 )(P / -P,; 



post 



Vn /pre r p r 

_P U — E B + ir\ 



(23) 



where ip is the three-dimensional TOPT analog of the Bethe-Salpeter wave function for the 
bound state, and R is regular at the bound-state pole. Inserting Eq. (^) into the second 
line of Eq. (|18|) , taking the residue at P° — > E^, and rearranging, we have 



^ B (Pl,P2,P3)=G SB ( Pl ,p 2 ,P3) / 



dp'i dp' 2 dp' 3 

(2tt) 3 (2vr) 3 (2tt) 3 



x 



E Vf B (p k , Pi, Pi! P*. P^ P,)(2vr) 3 5 (3) (p, - p}) 

3=1 



+ K) BB (pi,P2,P3;pi,P2,p 3 



^ B (P'1,P 2 ,P 3 ) 



where we have defined 

* Bb (Pi,P2,P3) = Mi(Pi)M2(p 2 )M 3 (P3) ^ b (Pi,P2,P3) 

and 



V; 



I rn k rrii 
tk{Pk)ei{Pi) 



^/ B (Pfc,P/,Pi;Pfc,Pl,Pi) Mfe(p' fc )uj(p{), spin-l/2 



Vf B (Pk,PhPj;p' k ,p'i,Pj 



£k(p' h h(Pi) 



spin-0. 



(24) 



(25) 



(26a) 



Vo 



«i(Pi)«2(P2)«3(P3) fKf B (Pi ; P2,P3;pi,P2,P 3 )] «i(p'i WP2WP3X spin- 1/2; 

V r BB (Pl,P2,P3;Pi,P / 2 ) P3 



ni\m2vn 3 



eilPije2lP2je 3 ^p3j 



ni\m2m 3 



ei(p / i)e2(p 2 )e3(P3) 



spin-0. 
(26b) 



Here the two-body interaction V^ B , which consists of a sum of two-particle irreducible time- 
ordered graphs between particles k and I {Vf 3 ), multiplied by spinor factors, depends upon 
the momentum of the noninteracting particle Pj. This is due to the term P° — ej(pj) in the 
denominator of the potential [see, for example, Eq. fllBDl. As noted earlier, the connected 
part vanishes by definition if there is a noninteracting particle. Therefore, cluster 
separability is unaffected by dependence of the two-body interactions that are internal to 
on the momentum of the spectator. The bound state for a two-body cluster in the limit 
that the third particle does not interact (i.e., for sufficiently short-range interactions and 
when there are no zero-energy bound states) derives from and the potentials internal 
to it have no dependence on the spectator momentum. 



Equation (24) is a Schrodinger-like relativistic equation: it is three dimensional, has a 
global relativistic propagator [Eq. (|T2"D], and reduces to the Schrodinger equation in the 
nonrelativistic limit. Also note that our interaction has energy dependence and hence retar- 
dation. 
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VI. CONCLUSION 



We have examined the full three-body Green's function for the case of one-boson- 
exchange interactions and positive-energy spin-0 or spin- 1/2 particles. We expanded out the 
Bethe-Salpeter equation for the Green's function into an infinite series of four-dimensional 
graphs. After performing a Fourier transformation of the internal energies into relative 
times, we integrated over the relative times, leaving an expansion of the full Green's func- 
tion which only involved three-dimensional internal variables (and hence integrations), while 
still depending upon the full four- dimensional nature of the initial and final states. 

Concentrating upon those graphs which contribute to the three-body bound-state pole, 
we re-summed the series into three factors: a three-dimensional Green's function G3 obeying 
an iterative equation, and pre- and post-factors which link the three-dimensional G3 to the 
off-shell states of the four-dimensional theory. The bound-state equation was then extracted 
from G3 and shown to have a Schrodinger-like structure involving a global relativistic prop- 
agator. The one-boson-exchange potential was shown to be augmented by factors which for 
spin-1/2 particles are plane- wave spinors, and for spin-0 particles are kinematical factors. 
Although the bound state is determined without reference to them, the pre- and post-factors 
are needed when interactions are considered. The current must include off-shell factors to 
account for the introduction of four-momentum into the graph through the interaction. The 
complete analysis of these currents, however, is left to future papers. 

Numerical calculations involving three bosons are under way. They are based upon 
Eq. (p4|) in the limit in which Vq — > 0, and compare the full (retarded) Vj to an instant 
approximation, as well as comparing these forms to those proposed by others. Full relativistic 
kinematics are used in conjunction with these relativistic interactions. 
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APPENDIX A: TIME- ORDERED PERTURBATION-THEORY RULES 

The rules for calculating graphs in a time-ordered perturbation theory using positive- 
energy particles are 

1. Assign an overall factor of 

(i) —i if no particles interact, 

(ii) 1 if only two particles interact, 

(iii) i if all three particles interact. 

2. Assign a factor of (2ti) 4: 5^[J2 (pj — Pi)}, where the sum is over interacting particles. 
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3. To each particle line with no interactions between the initial and final states, assign a 
factor of 

(2n)U^(p f -p t )N+(p t ) . 

4. To each final particle line emerging from it's last interaction, assign a factor of 

^ + (P/) _ 
p°j — e + ii] 

5. For each vertex on particle j, 

(i) Conserve 3-momentum, 

(ii) assign a factor of: 

ryv+( P ;) 
fa(\pj-p'j\) 

where p'j is associated with an earlier time than that with which pj is associated. 

6. For each unconstrained 3-momentum p, assign a factor of 

dp 



I 



(2nf 



7. To each time slice between vertices, assign a factor of: 

1 



P -E 1 -E 2 -E 3 -E m Um + iv ' 

where, 



E n = 



e » = yPn + m n-> ^ OT internal particles; 
pi, for external particles, 



and m ranges over the exchanged bosons existing during the time slice. 
8. For each initial line, a factor of 

1 



p® — e + 



12 



*3 

FIG. 3. An enlargement of the Feynman diagram on the left of Fig. ^. Here we also add the 
time and three-momentum labels used in expression (|B1|) to the particle lines. 



APPENDIX B: REDUCTION FROM FOUR TO THREE DIMENSIONS 

In Fig. ^| we have the Feynman diagram representing an interaction between particles 
one and two, "followed" by an interaction between particles two and three. As in Eqs. (|T^- 
|14D , we assume that this graph is embedded within other graphs. Using the notation of 
Fig. |3|, our positive-energy approximation, and performing a Fourier transformation of the 
time-component of momentum, we have 



dt f {e lP f l f {-i)6(t f - tje-^t-* 1 ) {-i)9{t f - t^e"^*'-^ {-i)9{t f - t^e^^) 
x J dU(-l)e- lP ? u J dt x dt 2 dt' 2 dt 3 
( iT 3 N 



2wi(p 3 - p' 3 ) 



(-i) 6(t' 2 - t 3 ) e -^(*2-*3) + 9 fe _ Q e +i"i(t> 2 -t 3 ) 



X I 



{-i)0{t' 2 - t 2 )N+e- ie '^~ t2) 



+ „-ie3(*3-ti) 



X 



2u 3 (p! - pi; 



-i)9(t 3 - tJN+e 

b 

(-i) 0(ti - t 2 ) e -^3(ti-t 2 ) + (-i) 0^~tij e +itJz{tl - t2) 



(Bl) 



x i 



-i)0(t! - U)N+e 



+ e -iei(tx—ti) 



{-i)9{t 2 -t t )N+e 



where the first two lines of expression ( |B1[ ) account for the fact that this graph is embedded 
within another graph, and we have labeled the exchanged-boson 9 functions A, B, a, and b. 
Using the identity 



(x - y)9(x' -y) = 9{x - x')9{x' -y) + 9{x' - x)9{x - y) 



(B2) 



one can expand the expression ( pi|) into the 12 terms shown in Fig. [|. The labels above the 
columns show from which 9 functions each column originates. The top four diagrams are 
three-particle reducible, and are therefore representable as iterations of the two-body force. 
The other eight diagrams are part of the three-dimensional three-body force. 
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Let us choose the third diagram in the (aA) column of Fig. |2|. This diagram results from 
the 9 functions 

0(t f - - t' 2 )9(t' 2 - h)9{h - t 2 )9(t 2 - U) . (B3) 

Defining the time interval variables 

T 5 = tf-tl, U = ti-t' 2 , T 3 =t' 2 -t 3 , T 2 = t 3 -t 2 , Ti=t 2 -U, T = ti, (B4) 

and noting that \d(tf, t X) t 2 , t 3 , t 2 , £j)/<9(r 5 , r 4 , r 3 , r 2 , r 1; r ) | = 1, Eq. ( [B~T| ) has the form (for 
this particular combination of 9 functions) 

dT 5 dndr 3 dr 2 dTidTo 9{t^)9{t a )9{t 3 )9{t 2 )9{t 1 )9{tq) 

/ r 2 r 3 x 



X 



,2^i(p 3 - p' 3 ) 

Ti r 2 



-io)i(r 3 ) 



-11^3 (T4+T3+T2) 



% 2W 3 (P1 -Pi)y 

Using the identities 

j°° dr ^ p h p ?) = 2n 5(P° - P\ 
J —00 



jy+ e _i<E i (t"4 +T3+T2 +n ) jy + e _ie 2 (n ) 



(B5) 



(B6a) 



dr 9{t) e 



iT(a~b+irf) 



a — b + ii] ' 

we can calculate the integrals in Eq. ( |B5| ). Simplifying the result, we have 

1 



ei - e 2 - e' 3 + £77 



1 



2cu 3 (p! - pi) P f - e i - 4 - e' 3 - ^3 + iv 

1 



x 



r 2 jv+ 



r 3 iv 3 + 



X 



2a; 1 (p 3 - p 3 ) Pf " ei - 4' - 4 " " ^3 + «7 A /2^(p 3 -p 3 ) 
1 r 2 iV+ 



P9 



J - e x - e 2 - e 3 - u 3 + vq 2u 3 ( Pl - p[ 



Pj -e 1 -e 2 -e 3 + ir) 



(B6b) 



(B7) 



This result agrees with the rules given in Appendix. [A[ The effect of the pre- and post- 
factors on this "internal" graph is to introduce the total energy for the initial and final 
states. Each time interval corresponds to a denominator with the on-shell energy of each 
existing particle subtracted from the total energy. The graph begins and ends with a Go 
factor. Note that the graph in Fig. (g) does not include these G factors. Instead it includes 
only those factors within the brackets. These factors are also the ones which contribute to 
V . 
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APPENDIX C: CLUSTER SEPARABILITY AND TOPT 



Cluster separability (CS) implies that if we describe particles propagating using a Green's 
function, and if one cluster of particles does not interact with another cluster of particles, 
then we can perform a separation of variables (between these two clusters) on the Green's 
function (Z§[2B|. 

Let us describe three distinguishable particles propagating with initial momenta pi,p 2 ,P3 
and final momenta pi,p 2 , p' 3 by the Green's function G(p' 1 ,P2,P 3 ]Pi,P2,P3)- CS tells us that 
in the absence of interactions, we must have 

G (pi > p'2 , p's ; pi > vi , pa) — ► £qi (pi ; pi ) G%2 (p'2 ; P2 ) ^; 3 } (ps ; Pa ) • ( c 1 ) 

If we have only two of the particles interacting (say 1 and 2), we have 

G{p , nP2iPziPiiP2,Pz) — ► Gi%(p'i,P2iPi>P2)QShP3iP3) ■ (C2) 
Finally, if all three particles interact, we have the fully connected Green's function 

G(Pl,P2,P3iPuP2iP3) * Gc ) (Pl,P2,P3iPl,P2,P3) • (C3) 

As these are the only possible cases, CS states that for the most general three-body Green's 
function we have Eq. flT6|) , which we restate here 

G (pi > p'2 , p's ; pi . P2 , P3 ) = G^i (pi ; Pi) G& (p'2 ; pi ) (Pa ; ps ) 

3 

+ Z! G%{p' P p' k ; Pj,Pk)G^ (p'i; Pi) 
1=1 

+ Gi 3) (p'l,P2,P3,Pl,P2,P3) • © 

We are working in momentum space, where the Green's function is simply the result of 
performing a Fourier transform on the configuration space Green's function. For the ra-body 
system, 

/n 
J] d 4 x t d% e^-^Gix,, r„ ://l //,. ). (C4) 
i=l 

Please note that some authors choose to note the Dirac 5 functions explicitly, by factoring 
them out of the momentum space Green's functions. In this case, the left hand side of 
Eq. (0) would be 

(2vr) 4 (5 (4) (pi + • • • + p' n -pi p n )g(p' 1 ,...,p' n ;p 1 ,...,p n ) = ■■■ . 

We have chosen not to perform this separation. 

1. The TOPT case 

Let us now examine the results of CS in the context of the three-body TOPT rules given 
in Appendix 0. 
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a. No interactions 



This case is the simplest, giving 



X 



Pi - e\ + it] 
(27r) 4 ^)(p 2 -^ 2 )iV+(p 2 ) 
P% ~ £2 + iV 



x (27r) 4 ^)(p 3 -^)iV+(p3) ^ 
P° - e 3 + if] 



which implies that 



gSKp[;pi) = igMiPi) = W^fa - P [) N+{pi ^. . (ce) 

Pi - ei + 



6. Two particles interacting 

For definiteness we will assume that particles 1 and 2 are interacting, and that particle 
3 is the spectator. This implies that of the four interactions V , V\, V 2 , and V3, only V 3 is 
nonzero. 

From the TOPT rules given in Appendix [A], we obtain the Green's function Q for the 
case of particles 1 and 2 interacting. It is possible to write it in operator notation (defined 
shortly) in the form 



G(Pl,P2,P3iPl,P2,P3) 



z(27r)V 4 )(p 3 -p 3 ) — 

p'i - e 3 + %r\ 



Pi -P2J 



Pi 



/o 



e[ + ir\ p' 2 



x 



vf D) iv 1 + iv 2 + 



+ 



e 2 + irj 

(3D^4D) 



1 



X 



Pi — ei + if] p 2 - e 2 + if] 

#S (P3i Ps) GSl2(Pl,P2,Pl,P2) 



(C7) 
(C8) 



where the last line identifies ^MPiiP^Pi^)- First we will define the different forms of 
the potential V3, and then we will define the global propagator for cluster I, G l , and the 
wave operator fl 3 . 

As shown in item [7| of the TOPT rules in Appendix |A], the denominator of each time 
interval takes a different form depending upon whether all, some, or none of the particles are 
external (e.g. initial or final) particles. We may factor the fully connected Green's function 
into parts based upon the forms of these denominators and their dependence upon the initial 
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and final relative energies. We use the superscript (4D) to denote the part which depends 
upon both the initial and the final relative energies, (4D — > 3D) and (3D — > 4D) to denote 
transitional parts depending only upon the final or initial relative energies, respectively, and 
(3D) to denote the part which is independent of both the initial and final relative energies. 
These forms are four-dimensional, transitional, and three-dimensional, respectively. For 
the case of these parts are simply different forms of the two-body interaction. (The 
parts are more complicated for the three body case, as we see later in this Appendix and in 
Sec. [IV].) We may see the relevant differences in these forms of the two-body interaction most 
easily by briefly restricting ourselves to the one-meson-in-flight approximation. Extension 
to the complete interaction is straightforward. 

The four-dimensional form of this interaction follows from both particles being external 
particles. Denoting the two-body center-of-mass energy as P® 2 = P\ + P®, 

V (4D) _ r\ i r . 



'MlPi - Pil) Pl ° 2 -Pi-Pz-u + iV J2u(\p 2 - p' 2 \ 



'Mil* - Si) p i°2 - p? - - u + *n \IM\pi - Pil) ' 



Note that the denominator depends upon either the initial or final energies of the particles 
p°, rather than the "on-shell" energy e. This is due to the fact that they are always in either 
initial or final states, and allows for the full four-dimensional nature of Q. This only appears 
in the Born term, in which the two particles interact only once. 

The transitional form of the potential has two variants: one initial particle with one 
internal particle, and one final particle with one internal particle. Here we show the second 
variant explicitly 



(4D^3D) 



ri i r 2 



2cu(|q - pil) F i2 - d - Pi -u + ir) y^QP - q - p' 2 |) 

+ Eh I Fi ( C10 ) 

/2w(|P - q - p' 2 |) p i°2 - Pi -e 2 -uj + ir) J2u(\q - p[\) ' 



Note that each denominator depends upon one particle's on-shell energy e, and the other 
particle's final energy p°. This is due to one particle going into its final state at the beginning 
of the interaction, while the other does not do so until the end of the interaction. The on- 
shell (internal) particle reflects our integration over the internal time variables, while the 
"off-shell" (final) particle reflects the (fully specified) four-dimensional nature of Q. The 
other variant of this form V^ 3D ^ 4D ^ is similar, and involves the initial particle states. 

Finally we have the three-dimensional form, where all of the particles are internal parti- 
cles, and hence we have integrated out their energy dependence. 

V r 3 (3D) - Tl 1 12 



2w(|q - q'|) P i2 -ei-e'2-u + ir] ^2u;(|q'-q|) 



+ . Fa -Q ^ . Fl (Cll) 

'2u(\q' - q|) P i2 - A - £2 - w + if] J2a;(|q - q'|) 
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Note that all of the particles are on shell, as none of them are initial or final particles. This 
reflects our integration over all internal time variables. This is the form which appears in the 
two-body bound-state equation, and is independent of the initial and final relative energies. 

We have only been considering the one-meson-in-flight approximation, but these com- 
ments hold for the general case of V3. 

We also define the two-body global propagator for cluster I 

c i - 1 
°" Pf k -e j -e k + i V ' 

where P° fc = p° + pi is the total energy of the pair. 

It is useful to define the right and left two-body wave operators for cluster j 

00 

it} = £ (v^NtNfGi) , (C12a) 

m=0 

00 

nf = £ (Givp N tK) m . (ci2b) 

m=0 

These wave operators transform the three-dimensional two-body free propagator G J into the 
full three-dimensional two-body Green's function G\: 

= Gi + GiV^N+N+Gi . (C13) 



c. All three particles interacting 

Finally, we have all three particles interacting. As in the two-body case [Eq. (|C7f) l, we 
may separate out the Born terms from the fully connected Greens function. In this context 
Born terms are defined as those in which all three-body reducible time intervals contain 
either an initial or final particle; there is no factor of Go, the three-dimensional three-body- 
reducible time interval. Examples of the three-body Born diagrams are given in Fig. |j. 
These graphs are the analog of the two-body V 3 ( ' 4D \ The simplest example is one of the 
terms associated with the graph in Fig. |](d); a single TOPT three-body force. Choosing 
the term analogous to that shown in the second diagram in column (bA) of Fig. ^|, but with 
external legs, we have 
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(a) 



(b) 



(c) 





(d) 



(f) 







(g) (h) 
FIG. 4. The TOPT graphs representing the 'Born' terms for the fully-connected three-body 

Green's function. The open circles represent one or more two-body forces, the boxes in graphs (b) 

and (c) represent two-body forces, and the boxes in graphs (d) through (h) represent three-body 

forces. The identifying feature is that no graphs contain a three-body reducible time interval that 

has no initial or final particles. 



- l (2,n^(p - n SM> jgte) j£W 

° Pi - e x (pi) + if] p' 2 ° - e 2 (p' 2 ) + ii) V 'i ~ e 3 (P3) + ^ 

r 2 iY 2 + (P-p' 1 -p 3 ) 1 



x- 



X 



2c^(p 3 - p 3 ) ^° - Pi - e 2 (P - p'i - Ps) - Ps ~ ^(P3 - Ps) + it] 
r 2 iV+(p 2 ) 1 r 3 iY 3 + (p 3 ) 



2cj( Pi - pi) P° ~ Pi -Pi- P'i - "(pi - P'i) - ^(Ps - Ps) + it] ^J2uj(p 3 - p 

i r 1 iv+( Pl ) 



po _ p f _ p o_ p o_ w ( pi _ p Q + ir] ^(pi - pi) 
1 1 1 



Pi - e i(Pi) +iV P§ - e a(P2) + itj pi - e 3 (P3) + it] 



(C14) 



Note that both initial and final particle energies are needed: this is a fully four- dimensional 
graph. 

In Sec. [V of the main text we show that when one excludes the three-body Born terms, 
can be factored in a similar way to Eq. (|C7| ), as shown in Eq. (|I7|), which we repeat 
here for convenience 

= i (2n)W\P f - P,)Cr 3D) Gf D) f§^ . (0) 
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Here the pre- and post-factors (/pr ( ?~ >4D ' ) and /p OS t^ ) correspond to the transitional factors 
containing V^ 30 ^ 40 ^ and V^ 4D ^ 3D ' ) in Eq. flC7|). They connect the three-dimensional G3 to 
the four dimensional external world. 

It is interesting to see how G)?> of Eq. fllTf ) vanishes in the limit that one particle does 
not interact (i.e., for sufficiently short-range interactions and when there are no zero-energy 
bound states). If particle 1, say, does not interact then V 2 , V 3 , and V vanish: by definition 
these potentials involve particle 1 interacting at least one time. This causes both / pre and 
/post [Eqs. (|19|)1 to vanish, as they both involve factors of two of the four possible interactions 
iy%i V2, V3, and Vo), three of which must now vanish. Note that this does not require G3 of 
Eq. (|18D to vanish, merely the factors which multiply it in the definition of Q^> . In fact, even 
in the absence of any interactions G3 has a nonzero value (Go), however it is not physically 
meaningful; if either f pre or f post vanishes, then the physically meaningful vanishes. 
Another way to state this is that is fully connected, and therefore must vanish in the 
limit that one particle does not interact. In the factorization of we have introduced, the 
pre- and post-factors ensure connectedness and therefore ensure that it vanishes. 
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